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We investigate the structural properties of nuclear star clusters in late-type spiral 
galaxies. More specifically, we fit analytical models to Hubble Space Telescope images 
of 39 nuclear clusters in order to determine their effective radii after correction for the 
instrumental point spread function. We use the results of this analysis to compare 
the luminosities and sizes of nuclear star clusters to those of other ellipsoidal stellar 
systems, in particular the Milky Way globular clusters. Our nuclear clusters have a 
Q\ ■ median effective radius of r\T = 3.5 pc, with 50% of the sample falling between 2.4 pc < 

r e < 5.0 pc. This narrow size distribution is statistically indistinguishable from that of 
. Galactic globular clusters, even though the nuclear clusters are on average 4 magnitudes 

brighter than the old globulars. We discuss some possible interpretations of this result. 
Q-f From a comparison of nuclear cluster luminosities with various properties of their host 

galaxies, we confirm that more luminous galaxies harbor more luminous nuclear clusters. 
It remains unclear whether this correlation mainly reflects the influence of galaxy size, 
mass, and/or star formation rate. Since the brighter galaxies in our sample typically 
have stellar disks with a higher central surface brightness, nuclear cluster luminosity 
also correlates with this property of their hosts. On the other hand, we find no evidence 
for a correlation between the presence of a nuclear star cluster and the presence of a 
large-scale stellar bar. 

Subject headings: galaxies: nuclei — galaxies: spiral — galaxies: star clusters 
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1. Introduction 

In most formation scenarios for spiral galaxies, the central bulge is the "trashbin of violent 
relaxation" where a dynamically hot stellar component has formed either through external potential 
perturbations such as early mergers of proto-galaxies (e.g. Carlberg 1992), or perhaps via internal 
effects such as violent bar instabilities (Norman et al. 1996). While this "nature or nurture" debate 
is far from settled, the recent discovery of the so-called M-a relation (Gebhardt et al. 2000; Ferrarese 
& Merritt 2000) has made it clear that the structure of the very nucleus, i.e. the black hole and its 
immediate vicinity (<S lpc), is closely linked to the overall dynamical structure of the stellar bulge. 
The latest-type spirals, then, must have lived very sheltered and uneventful lives, since their central 
"trashbin" is virtually empty. Indeed, these disk-dominated, bulge-less galaxies (e.g. Matthews &: 
Gallagher 1997) often have slowly rising rotation curves that indicate a nearly homogeneous mass 
distribution on scales < 1 kpc. On these scales, gravity therefore hardly provides a vector pointing 
at the center, and it is not obvious that the nucleus of these galaxies is well-defined and a unique 
environment. 

It therefore came as a surprise that recent observations with the Hubble Space Telescope (HST) 
have shown that the photocenter of many late- type spiral galaxies is occupied by a compact, 
luminous stellar cluster (Phillips et al. 1996; Carollo et al. 1998; Matthews et al. 1999). In a recent 
paper (Boker et al. 2002, hereafter Paper I), we have shown that at least 75% of spiral galaxies 
with Hubble types later than Sc harbor such a nuclear star cluster. The fraction of nuclear clusters 
(NCs) appears to be somewhat lower in earlier Hubble types, as indicated by the results of Carollo 
et al. (1998). However, smaller sample sizes and the observational difficulties due to the more 
complex nuclear morphologies of early-type spirals and the decreasing contrast between a central 
star cluster and a luminous bulges make this result somewhat uncertain. Typical luminosities of 
NCs are in the range 10 6 - 10 7 L (Paper I) which makes them much brighter than average stellar 
clusters in the disks of nearby spiral galaxies (Larsen 2002), and comparable to young "super star 
clusters" in luminous merger pairs (e.g. Whitmore et al. 1999) or circumnuclear starforming rings 
in spiral galaxies (e.g. Maoz et al. 2001). 

Despite the recent progress, the formation mechanism of NCs remains largely a mystery. So 
far, no satisfying explanation has been put forward to explain the high gas densities that must have 
been present to enable the in-situ formation of such compact and luminous objects in the nuclei 
of these disk galaxies with shallow potentials. Plausible alternatives to in-situ formation include 
the infall of evolved clusters, possibly combined with more recent star formation events due to the 
accretion of gas in the galaxy nucleus. In order to gather evidence for or against any of these 
scenarios, it is essential to identify the stellar population(s) that comprise the NC, and thus to 
constrain the cluster formation history. Another crucial piece of information in this undertaking is 
an accurate estimate for the cluster mass. 

The stellar masses of NCs can be derived - under the assumption of spherical symmetry - if one 
has accurate knowledge of (a) the stellar velocity dispersion and (b) the cluster light distribution. 
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This kind of analysis has been successfully applied to NCs (Boker et al. 1999) as well as young 
clusters in the disks of nearby galaxies (Smith & Gallagher 2001; Mengel et al. 2002). Our team 
has obtained high-resolution spectroscopy of ~ 15 clusters from the sample of Paper I. These 
data will address point (a) above, i.e. to accurately measure the stellar velocity dispersion of the 
clusters (Walcher et al. 2003, Walcher et al. 2004, in prep.). The purpose of the present paper is 
to discuss point (b), namely the size and spatial structure of the NC sample of Paper I. Because 
these objects are only barely resolved even with HST, one has to carefully correct for the effects of 
the instrumental point-spread function (PSF). Because standard deconvolution algorithms are not 
well suited in the realm of "almost point sources", we have used the specialized software package 
ISHAPE (Larsen 1999) to find the best match between the data and PSF-convolved analytical 
models for the cluster shapes. We describe the details of this exercise in § 2, and compare the 
resulting structural parameters of NCs to those of other cluster populations in § 3. 

Because nuclear star clusters are a relatively new aspect of galaxy morphology, it seems a 
logical approach to search for correlations between the properties of NCs and those of their host 
galaxies. As has been shown by Carollo (1999) for spirals of earlier Hubble type, such empirical 
correlations might reveal clues about the mechanism(s) that regulate the gas supply in the very 
centers of spiral galaxies, and will certainly improve our understanding of nuclear activity in spirals. 
Moreover, because the presence of a NC appears to be a common feature in all types of spiral 
galaxies, the dynamical processes that lead to nuclear starbursts will be an important contribution 
to the discussion on bulge formation scenarios and a possible evolutionary connection between the 
various Hubble types. For these reasons, we present in § 4 a number of correlation plots between 
NC sizes/luminosities and various host galaxy properties such as size, Hubble type, HI mass, star 
formation rate, or total magnitude. We summarize our results in § 5. 



2. Structural Analysis of Nuclear Star Clusters 
2.1. Methodology 

Star clusters in all but the closest galaxies are marginally resolved sources, even with the 
excellent resolution of HST images. In order to measure reliably their intrinsic sizes and structural 
parameters, it is therefore crucial to have both observations with high signal-to-noise ratio (SNR) 
and accurate knowledge of the instrumental point-spread function (PSF). Insufficient sampling of 
the PSF can complicate matters further. For example, the pixel size of the WFPC2 (0.1" and 
0.0455" for the WF and PC chips, respectively) in general does not provide Nyquist sampling 
of the HST resolution element R = X/D. While sub-pixel dithering can alleviate this problem 
to some extent, HST observations of extragalactic star clusters in all but the closest galaxies are 
almost always confined to the regime of "marginally resolved" . 

In this situation, none of the conventional deconvolution methods such as the Maximum En- 
tropy Principle (Burch, Gull, & Skilling 1983) or the Richardson-Lucy algorithm (Richardson 1972; 
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Lucy 1974) work reliably. Instead, differential aperture photometry (e.g. the magnitude difference 
mi — 7714 between concentric apertures of 2 and 4 pixels, Holtzman et al. 1992) of stellar clusters 
and comparison to individual stars broadened by Gaussians has been used to estimate cluster sizes. 
However, a Gaussian is generally a poor approximation to the surface brightness profile of stellar 
clusters. In particular, the wings of a Gaussian fall off too rapidly, and cannot adequately describe 
the outer regions of either old globular clusters (GCs) in the Milky Way or young GCs in external 
galaxies. 

More recently, a number of authors (e.g. Kundu & Whitmore 1998; Larsen 1999; Carlson k, 
Holtzman 2001) have independently developed techniques to determine more reliably the structure 
of stellar clusters in external galaxies. The various methods are conceptually similar in that they 
assume a parametric model for the intrinsic light distribution. A particular model is considered to 
be a good description of reality if - after convolution with the instrumental PSF - it matches the 
observations. The "best" model is the one that minimizes the x 2 -difference to the data. The most 
widely used analytic models for the projected density profile of stellar clusters are those of King 



Here, z is defined by z 2 = a\x 2 + aiy 2 + a^xy, with x and y denoting the coordinates relative to the 
center of the profile and constants cti, 02, and CI3 which depend on the major axis, ellipticity and 
orientation of the model. The so-called concentration index c is the ratio between the tidal radius 
r t and the core radius r c . 

In the regime of marginal resolution, the results of Carlson & Holtzman (2001) have shown that 
unless the SNR of the observations is exceptionally high (> 500), it is very difficult to distinguish 
between King models of different concentration indices. On the other hand, they find that the 
effective radius r e , i.e. the radius which contains half the cluster light in projection, can be fairly 
robustly estimated. The value of r e does not depend strongly on the exact shape of the assumed 
model. This is in agreement with the results of Kundu & Whitmore (1998) and has been extensively 
tested by Larsen (1999). In what follows, we will describe our attempts to measure the effective 
radii of the NCs found in the survey of Paper I. 

For our analysis, we have made use of Soren Larsen's software package ISHAPE, which is 
described in detail in Larsen (1999). This software package has been well tested in a number of 
studies (e.g. Mengel et al. 2002; Larsen et al. 2002b), and has proven to be a robust and user- 
friendly tool for structural studies of compact star clusters. In brief, the program convolves a two- 
dimensional analytical model (e.g. Gauss-, Hubble-, Moffat-, or King-profiles) with a user-provided 
PSF, and finds the particular parameter set that best describes the observations, including the 
exact location of the cluster centroid on sub-pixel scales. Comparison with the pure PSF (i.e., a 
PSF-convolved <5-function) allows to assess whether the clusters are resolved at HST resolution. 
Because r e is not defined for functions with divergent integrals such as Hubble-type profiles, we 
have only attempted to fit the data with the above-mentioned King models as well as Moffat profiles 
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which are defined as 

/W = (TTW • < 2 > 

with z defined as for the King model, and i > 1 denoting the so-called power index. 

In order to construct accurate PSFs for our analysis, we have used the TinyTim software 
(Krist k, Hook 2001), which accounts for the spectral passband used, the plate scale variations 
across the field of view, and the spectral energy distribution (SED) of the source. Although all 
of our images were taken through the same filter (F814W), we have constructed a separate PSF, 
according to the exact position of the NC on the detector array. We note, however, that this only 
marginally affects the results because the WFPC2 PSF varies only slightly across the PC field. 
The choice of the correct SED for constructing the PSF is somewhat more critical in estimating 
the NC sizes. Unfortunately, spectroscopic information is available for only a limited number of 
NCs. Early results of our spectroscopic follow-up survey (Walcher et al. 2003) show that the light 
of many NCs is dominated by a stellar population in the age range 50-300 Myr. This agrees with a 
number of case studies of individual clusters, e.g. in M 33 (Gordon et al. 1999), NGC 2403 (Davidge 
& Courteau 2002), IC342 (B6ker et al. 1999), or NGC 4449 (B6ker et al. 2001). We therefore chose 
the SED of an A5 star for constructing the TinyTim PSF for all objects. We have verified that 
an alternate choice for the cluster SED, e.g. that of a cool giant with spectral type M3, does not 
significantly change the ISHAPE results. 



2.2. Results 

We have used ISHAPE to fit all 59 NCs of Paper I within the aperture radius that is dominated 
by emission from the NC (at least 5 r c ). Some typical examples that illustrate the ISHA PE results 
are presented in Figure 1. We used only circularly symmetric models within ISHAPE because 
visual inspection of the fit results confirmed that the majority of clusters show little evidence for 
flattening. In a few cases, the clusters are asymmetric or have a complex morphology, most likely 
due to dust lanes or faint companion clusters. These objects were excluded from the ISHAPE 
analysis. In addition, we excluded those clusters with insufficient SNR to obtain reliable fits, i.e. 
with SNR < 40 (Larsen 1999). 

The remaining 39 clusters were fit using King profiles with c = r t /r c values of 15, 30, and 
100, as well as Moffat profiles with power indices of i = 1.5 and i = 2.5. In principle, the galaxy 
light underlying the NC must be subtracted before fitting the cluster profile. The ISHAPE code 
simply subtracts a constant background level, measured at the egde of the fitting aperture. While 
this algorithm certainly is less than perfect, and more sophisticated methods to fit the surface 
brightness profile of the galaxy disk could be employed, we have verified that the resulting values 
for r c are not systematically biased by this shortcoming. In fact, the cluster magnitudes calculated 
by integrating the best-fit ISHAPE model agree well with those listed in Paper I which were derived 
with a more careful estimation of the galaxy background (see § 3.1). The median (mean) difference 
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between these two sets of numbers is 0.07 (0.11) magnitudes when integrating the models over the 
same aperture used for the Paper I photometry. 

The best-fitting models and the resulting effective radii, both in angular and physical dimen- 
sions, are summarized in Table 1. The majority (29/39) of NCs have r e < 5pc, with median (mean) 
values of 3.5 pc (5.1 pc). As discussed further in § 3, these values are comparable to the typical size 
of old globular clusters in the Milky Way and other nearby galaxies. We emphasize that all clusters 
are resolved, i.e. models that simply rescale and position the TinyTim PSFs are inconsistent with 
the data by a large margin. For completeness, we also list in Table 1 the best-fitting analytic 
model. We caution, however, that in some cases, other models may also be consistent with the 
data. However, this has little impact on the uncertainties of the r e measurements, as explained 
further below. Our goal is not to determine the detailed shape of the cluster profile, but merely to 
obtain a robust estimate for its size. 

The uncertainties listed in Column 7 of Table 1 represent the range of r c for the best-fitting 
model profile that yields acceptable fits to the data. Here, "acceptable" is defined by the 3a 
confidence limit for one free parameter, i.e. Ax 2 < 9 (e.g. Press et al. 1992). Our analysis is 
conservative in the following sense: before evaluating the quality of the model fits, we rescaled all 
X 2 values such that the best ISHAPE model for a given galaxy had a x 2 value equal to the number 
of degrees of freedom, i.e. the number of data points (pixels in the aperture) minus the number 
of free parameters. This rescaling was necessary because ISHAPE appears to underestimate the 
intrinsic errors in the cluster images. As a result, many of the best x 2 values do not formally 
represent "good" models according to x 2 -statistics. As the scale factor was always smaller than 
one, the x 2 rescaling is equivalent to increasing the noise level of each data point or pixel, and 
consequently broadens the range of models that appear consistent with the data. Unfortunately, 
this rescaling of errors made a more rigorous statistical comparison between the different analytic 
profiles impossible. 

To illustrate the analysis, Fig. 2 shows for the example of NGC 2805 the rescaled x 2 values 
as a function of r c for all models considered. In this case, a Moffat profile with a power index of 
1.5 (MOFFAT15) and effective radius r c = 0.056" ± 0.003" yields the best fit to the data. The 3a 
confidence limit and the corresponding range of r c is indicated by the horizontal line in the inlay 
panel. As mentioned before, the effective radius is a rather robust quantity, i.e. different models 
yield similar values for r e , even if they are not formally consistent with the data according to the 
above criterion for Ax 2 . In the case of NGC 2805, this can be seen in Fig. 2: the "best" values 
for r c of all five models are within 20% of each other which is also the typical uncertainty listed in 
Table 1. 
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3. Discussion of Structural Parameters 

3.1. Comparison to Milky Way Globular Clusters 

A natural question to address with the data we have in hand is that of a possible size-luminosity 
relation for NCs. Since NCs in late-type disks are potential candidates for the seeds of bulge 
formation, they might follow the mass-radius correlation established for more massive ellipsoidal 
stellar populations, i.e. elliptical galaxies and bulges of intermediate- and early-type spirals (e.g. 
Kormendy 1985; Bender et al. 1992). On the other hand, the old GCs in the Milky Way - which 
constitute the best-studied family of stellar clusters - are remarkably free of any such trend: the 
effective radii of Galactic GCs are independent of cluster luminosity or mass (van den Bergh et 
al. 1991; Djorgovski & Meylan 1994; McLaughlin 2000). The same is true for GCs in the nearby 
galaxies M31 (Barmby et al. 2002), M33 (Larsen et al. 2002a), and NGC5128 (Harris et al. 2002). 

In Fig. 3, we plot the absolute I-band magnitude of the NCs, Mf , against their effective radius, 
r e . The procedure for measuring Mp is described in detail in Paper I. In brief, the challenge is to 
determine the light contribution of the stellar disk underlying the NCs. For this, we use two limiting 
models for the inward extrapolation of the galaxy disk profile. The Mf values in Column 5 of 1 
represent the average of the resulting two numbers; half their difference determines the size of the 
error bars in Fig. 3. Aside from a handful of bright clusters with effective radii larger than average, 
the bulk of the sample — 28 of 39 clusters — is found in the narrow range 1 pc < r e < 5 pc, while the 
clusters range over a factor of ~ 100 in luminosity, from —9 > Mi > —14. As we discuss below, 
relating luminosity to mass is nontrivial for NCs, but taken at face value, the situation illustrated 
in Fig. 3 is more reminiscent of GCs than of the bulges in early-type spirals. 

It is not only the absence of a correlation between size and luminosity that appears to connect 
the NCs in our sample to GCs. As noted in § 2.2, the distribution of r e for our NC sample is 
quantitatively similar to that of Galactic GCs. This is shown explicitly in Fig. 4, where we also 
compare the distributions of I-band magnitude for our NC sample and for 111 Milky Way GCs 
without an obvious or suspected core-collapse morphology, i.e., with well-determined r c values. 7 
Evidently, although the median Mj of NCs is fully 3-4 mags brighter than that of the Milky Way 
GCs, the widths of the two luminosity distributions (which corresponds to a distribution in mass 
for the GCs) are comparable. The median effective radius of the NCs is r~ e = 3.5 pc (interquartile 
range 2.4-5.0 pc) and that of the GCs is Yl = 3.3 pc (interquartile range 2.4-5.2 pc). A two-sample 
Kolmogorov-Smirnov test confirms statistically the visual impression that the two r e distributions 
are identical, with a probability Pks = 0.75 that they have been drawn from the same parent 



7 The Milky Way GC data are taken from the February 2003 version of the online catalogue compiled by Harris 
(1996), which contains 141 clusters with known values for angular half-light radius and heliocentric distance, King 
concentration index, absolute V-band magnitude, and (B — V) color and excess. For 96 of these, Harris also gives an 
apparent (V — I) color, which we correct according to E(V — I) = 1.25 E(B — V) (Cardelli et al. 1989) in order to 
derive Mi from My . We simply assign the average (V — I)o = 0.95 to the other 45 clusters. 
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distribution. We further note that the GCs in M31, M33, and NGC 5128 have essentially the same 
median r e as the clusters in Fig. 4, and the same independence from cluster luminosity (see the 
references above). 

Figure 5 provides an alternate illustration of this basic point, in a representation that allows 
for direct comparison with larger ellipsoidal stellar systems and comes closer to the "fundamental 
plane" typically plotted for such systems. Here, we show the average surface brightness, I c = 
Li/(27rr-g), vs. effective radius r c for the NCs, the 111 Galactic GCs from Harris' (1996) catalogue, 
and a number of E galaxies, spiral bulges, and dwarf Ellipticals (dE's) from the compilation of 
Burstein et al. (1997). The position of Milky Way GCs and NCs to the left of the plot shows 
again that the latter have essentially the same size as the former, but are substantially brighter 
and therefore at higher surface brightness. 

Of additional interest is the disconnect between the GCs and NCs on the one hand and the 
elliptical galaxies and spiral bulges on the other, a dichotomy first noted by Kormendy (1985). As 
we mentioned at the beginning of this discussion, the bulges and galaxies show a size-luminosity 
relation, while GCs do not. Equivalently, GCs have an enormous spread in surface brightness, while 
bulges and ellipticals at a given r e appear much more uniform. Although this fact is still not well 
understood, the alignment of our NCs with the GCs in Figure 5 would seem to underline a closer 
link between NCs and GCs, as opposed to one between NCs and classical galaxy bulges. 

For completeness, we have included in Figure 5 data for the nuclei of 5 nucleated dE galaxies 
in Virgo, from the sample of Geha et al. (2002). As those authors have already demonstrated, with 
a different version of the "fundamental plane", the dwarf-elliptical nuclei are also more similar to 
star clusters than to galaxy bulges. 

3.2. Possible Interpretations 

The apparent similarity between NCs and GCs could be an important clue to the nature and 
origin of the NCs themselves. First, however, their identification as an extension of the GC sequence 
to brighter magnitudes (Fig. 4) and higher surface brightness (Fig. 5) needs to be understood in 
terms of mass rather than luminosity alone. Of particular concern in this context is the interpre- 
tation of the spread of NC absolute magnitudes in Fig. 3 and Fig. 4. We consider two extreme 
scenarios — neither of which we mean necessarily to advocate as plausible per se, but which may 
serve to delimit further discussion: 

(1) If the stellar populations of all NCs shared a common age, then the five- magnitude spread 
in Mf would translate directly to a range of a factor of ^100 in NC mass, intriguingly similar 
to the corresponding variation among old GCs (Fig. 4). Were Fig. 5 then to be re-plotted using 
mass surface densities, the NC locus would simply shift vertically, by an amount depending on the 
typical NC age (through the /-band mass-to- light ratio). This would strengthen the connection of 
NCs to GCs and sharpen their distinction from galaxy bulges. 
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(2) If instead all NCs had roughly the same stellar mass, then the spread in Mf would have 
to arise from a range in luminosity-weighted age, with fainter cluster magnitudes (corresponding to 
systematically older ages) being purely the result of stronger population fading. If this were true, 
then Fig. 5 re-drawn in terms of mass surface densities would have the NC locus contracted into a 
relatively narrow line of slope —2. The small variation in NC effective radius would then be largely 
the result of a small mass range in our sample, and it would reveal little about the (non) existence 
of a mass-radius relation among these clusters. 

In order to gain some quantitative feel for the viability of these two alternatives, we have used 
the PEGASE population-synthesis code (Fioc & Rocca-Volmerange 1997) to calculate the evolution 
of /-band mass-to-light ratio and absolute magnitude for a stellar cluster of mass 7x 10 5 Mq, formed 
in a single burst of star formation with the IMF of Kroupa et al. (1993) between lower- and upper- 
mass cut-offs of 0.1 Mq and 120 Mq. The results, for a range of assumed metallicities, are presented 
in Fig. 6. Note that the mass-to-light ratio in the top panel is independent of the assumed total 
mass of the cluster, while the /-band magnitude in the bottom panel is of course sensitive to this 
parameter: for a given metallicity, the absolute magnitude of a single-burst stellar population of 
any mass m is obtained, as a function of age, by adding —2.5 log(m/7 x 10 5 Mq) to the appropriate 
curve in Fig. 6. The calculations for Yj and Mf both take into account the mass lost from a cluster 
due to stellar evolution, i.e., winds and supernova explosions. 

If the NCs were in fact a roughly coeval population of clusters with a narrow range in mass-to- 
light ratio Tj = M/Lj (scenario (1) above), then they could in principle have masses similar to the 
Galactic GCs if Tj were some 15 times smaller for NCs than for GCs (to explain the 3-magnitude 
difference in the peaks of the Mp distributions of Fig. 4). The top panel of Fig. 6 shows that a 
"typical" GC, with an age of 13 Gyr and [Fe/H] = -1.5, has Y 7 ~ 1.75 M L" 1 in these models. 8 
The same graph then shows that the NCs would have to be of order ~ 10 7 years old (to within 
a factor of two or so) in order to have Tj ~ 0.1 Mq Lq 1 , and thus a mass distribution similar to 
that of the GCs. Interestingly, this would make them potential analogues of the super star clusters 
forming in many nearby interacting and starburst galaxies (e.g., Whitmore et al. 1999). In this 
context, it is worth noting that these super star clusters tend to have effective radii r e ~ 2 — 10 
pc that do not appear to depend strongly on luminosity — reminiscent of NCs and GCs both (e.g., 
Barth et al. 1995; Zepf et al. 1999). Nevertheless, it is difficult to imagine how such a large fraction 
of late-type disks could have conspired simultaneously to form nuclear clusters such a short time 
ago. 

If, on the other hand, scenario (2) were correct, and the spread in NC luminosity reflected an 
age spread among single-burst stellar populations with essentially a single mass, then the bottom 
panel of Fig. 6 suggests that this common mass would have to be near 7 x 10 5 Mq, and the age range 
exceedingly broad: from 10 7 -10 10 years to just barely reproduce the observed — 14 < Mf < —9 



8 The observed number is about 1.2 MqLq 1 , which follows from an average F-band value of 1.45 in solar units 
(McLaughlin 2000), and a typical V -I color of 0.95 (see § 3.1). 
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(Fig. 4). (A change from the "fiducial" mass in Fig. 6 by a factor of even 2.5 in either direction 
would shift all the Mf curves up or down by a full magnitude, leaving an incomplete overlap with 
the observed NC magnitude distribution.) 

Even as oversimplified as these two scenarios are, it is not possible to discriminate further 
between them without additional information. We have therefore obtained short-wavelength spectra 
for a number of the NCs in the current sample, which we will fit with population-synthesis models 
to constrain the clusters' ages and star- formation histories, and from which we will extract stellar 
velocity dispersions to constrain dynamical mass-to- light ratios. Although this analysis is still in 
progress, the first results (Walcher et al. 2003) suggest a picture that is, not surprisingly, rather 
more complicated than either of the extremes just discussed. 

Walcher et al. (2003) show that the light-dominating population of stars in 10 relatively bright 
NCs range in age from ~ 10 7 -10 9 years old. That is, there is in fact a spread in (luminosity- 
weighted) age among the NCs — although it does not appear sufficient to explain the full range of 

given the single-burst models in Fig. 6, unless there is also a significant (order-of-magnitude) 
spread in the cluster masses. More importantly, however, Walcher et al. (2003) show that the 
single-burst hypothesis itself is almost certainly incorrect: the spectra of their sample of clusters 
show some direct evidence for a mix of stellar ages, and the dynamical mass-to-light ratios they 
derive are consistently higher than the value of T/ inferred from population-synthesis modeling, 
suggesting that an underlying old (faint but massive) population of stars may be present. 

In view of this, and given the similarity between GC and NC effective radii, an obvious first- 
order enhancement over the single-burst idea is the hypothesis that NC masses are in fact dominated 
by old stellar populations - perhaps bona fide GCs that have been dragged by dynamical friction 
to the centers of their host galaxies - but that their luminosities, and hence the young ages implied 
by their spectra, are due mainly to more recent star formation triggered by the infall of a smaller 
mass of circumnuclear gas. A more thorough exploration of this scenario is beyond the scope of 
this paper and must await the full analysis of our spectroscopic data (Walcher et al. 2004, in 
preparation). Nevertheless, it is worth noting that recent high-resolution interferometric studies of 
the nuclear gas dynamics in nearby galaxies have shown that molecular gas flows can indeed reach 
to within the central few pc of the nucleus (e.g. Schinnerer, Boker, & Meier 2003) and thus provide 
a possible mechanism for the fueling of such nuclear starburst events. 



4. Correlations with Host Galaxy Properties 

4.1. Luminosity and Mass of the Host Galaxy 

The work of Carollo et al. (1998) has shown that the luminosity of nuclear star clusters in early- 
and intermediate type spirals correlates with the total luminosity of the host galaxy. A similar 
result has been established by Lotz et al. (2001) for dE galaxies. In Figure 7a, we demonstrate 
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that such a correlation is also found in spirals at the late end of the Hubble sequence. Here, we 
plot the apparent I-band magnitude of the NCs, mf (from Paper I), against the total apparent 
blue magnitude, m^, of the host galaxy (Paturel, Bottinelli, & Gouguenheim 1994) as listed in the 
Lyon-Meudon Extragalactic Database 9 (LEDA). The relation still holds, albeit with a somewhat 
larger scatter, when plotting the absolute magnitudes Mf and Mg (Fig. 7b). The conversion 
from apparent to absolute magnitudes and the correction for Galactic reddening was performed 
using the foreground extinction and galaxy distances listed in Table 1 of Paper I. The solid lines 
in both panels indicate the best (least-squares) linear fit to the datapoints. The fit coefficients 
as well as their statistical significance (as measured by the Spearman rank-order coefficient) are 
listed in Table 2. The fact that despite the large intrinsic variations of the cluster luminosities 
discussed in §3.2, there still is a rather good correlation with galaxy luminosity suggests that the 
global properties of the host galaxy play a significant role for the formation processes of NCs. 

One possible mechanism that could link the occurrence of nuclear starbursts to the total galaxy 
luminosity are tidal interactions or minor mergers, because these can both enhance the global star 
formation rate (SFR) and cause significant infall of matter towards the nucleus via re-distribution 
of angular momentum. However, as demonstrated in Fig. 8, NC luminosity appears only weakly 
(if at all) correlated with the IRAS far- infrared (FIR) magnitude Mfir 10 which is often used as a 
proxy for the SFR in galaxies (see e.g. Telesco 1988, for a review). Assuming that tidal interactions 
contribute significantly to the SFR and hence the FIR luminosity in late- type spirals, this lack 
of a strong correlation between SFR and NC luminosity suggests that tidal interactions are not 
the driving force for the evolution of the very nucleus. This is supported by the WFPC2 images 
of our galaxy sample (Paper I) which show little evidence for recent star formation activity in 
their (central) disks. Instead, the smooth, shallow, and undisturbed stellar disks of many late-type 
spirals indicate a rather uneventful dynamical history, yet most of them have a luminous NC. 

An alternative explanation for the strong correlation between Mf and Mg could be that 
the properties of the NC are more governed by galaxy mass than by the global star formation 
rate. In principle, HI linewidths could be used to address this possibility: the inclination-corrected 
rotational component of the HI linewidth 

W^ = (W 20 -W rand )/smi , (3) 

based on the correction for random motions of Bottinelli et al. (1983), can be used as a proxy 
for the dynamical mass of the galaxy. Unfortunately, our galaxy sample was selected to have 
inclinations close to face-on and consequently the sin i corrections are highly uncertain. For the 
sake of completeness, we have netherless calculated the galaxy masses with the above expression. 



9 http: / /leda.univ-lyonl .fr / 

10 The apparent far-infrared magnitude, taken from LEDA, is calculated according to mpiR = — 2.51og(2.58 x f 6 o + 
fioo) + 14.75 where f6o and fioo are the IRAS fluxes at 60 /im and 100 /im (in Jansky). This expression is identical to 
the one used in the RC3. The distances listed in Paper I are then used to calculate M F ir. 
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Following Matthews, van Driel & Gallagher (1998), we adopted Wrand = 20kms~ 1 for our galaxies, 
which corresponds to a line-of-sight velocity dispersion a z = 5.5 kms -1 (Rhee 1996). However, 
using the inclination i as listed in LEDA, we found that for our sample, the HI width does not 
even correlate with total galaxy luminosity which is a contradiction of the Tully-Fisher relation. 
We thus conclude that for our sample of face-on spirals, is not a reliable indicator of galaxy 
mass, and that we therefore cannot draw any meaningful conclusions about possible correlations of 
NC properties and host galaxy mass. 

4.2. Central Surface Brightness of the Galaxy Disk 

Recent analysis of large galaxy samples (e.g. Blanton et al. 2002) has put on a solid statistical 
ground the notion that there are strong correlations between a number of physical properties of 
galaxies. In particular, Blanton et al. (2002) have shown that for disk-dominated, "exponential" 
galaxies (i.e., with Sersic (1968) shape parameter n < 1.5), surface brightness correlates with total 
luminosity. We confirm this finding for our sample, as demonstrated in Figure 9a. Here, we have 
plotted Mg, the total blue magnitude of the host galaxy, versus [if, the central I-band surface 
brightness of the stellar disk underlying the NC. The values for are listed in Column 3 of 
Table 1; they were measured from our WFPC2 images by averaging the two inward extrapolations 
of the surface brightness profile shown in Fig. 3 of Paper I at a radius of 0.024" (i.e. the border of 
the central WFPC2 pixel). 

Given that NC luminosity Mf correlates with host galaxy luminosity (Fig. 7b), and that the 
latter correlates with /Jj, it is not surprising that Mf is also correlated with /J}, as demonstrated in 
Figure 9b. Again, the fit coefficients and significance levels of both correlations are listed in Table 2. 
While it is true in general that higher surface brightness disks make it more difficult to detect faint 
clusters, we emphasize that this effect does not result in a significant selection bias. In fact, most 
galaxies without clear identification of a NC have /Jj values below average (Figure 4 of Paper I). 
Note that the values of fi® are as observed, i.e. they have not been corrected for inclination. 

So far, we have shown that more luminous galaxies have higher surface brightness disks as 
well as more luminous NCs. Is this simply a reflection of the popular notion that "bigger galaxies 
have more of everything"? While this is certainly part of the answer, it is likely not the whole 
story. This is demonstrated in Fig. 10, which plots Mf , as well as the cluster-to-disk luminosity 
ratio C/D, as functions of the physical size of the galaxy disk, i.e., the major axis diameter of the 
/xb = 25 isophote as listed in LEDA. (The disk luminosity is calculated as total galaxy luminosity 
minus NC luminosity, but with C/D < 0.01 in general the disk light is not significantly different 
from the total galaxy luminosity.) 

Apparently, Mf depends weakly on galaxy size: the Spearman rank-order correlation coeffi- 
cient in Fig. 10a is s = 0.36, with slightly less than 3-a significance. However, this actually reflects 
the stronger and more significant correlation between NC luminosity and galaxy luminosity (§4.1), 
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coupled with a similarly significant dependence of disk size on total disk luminosity: i?disk ~ -^disk 
in this sample. From Table 2 we have that NC luminosity scales with galaxy (disk) luminosity as 
L c \ ~ -^disk) anc ^ thus we can expect C/D = (A;i/Alisk) ~ -^disk 22 ~ ^disk 69 even if there is no direct 
connection between NC luminosity and global disk size. This scaling is indeed consistent with the 
weak trend in Fig. 10b. 

To first order, the NC luminosity is a measure of the efficiency with which gas is funneled 
towards the nucleus and converted into stars, regardless of the exact star formation history of the 
NC. This is true even if the cluster luminosity is dominated by a minor (in mass) but recent burst, 
because on average, younger populations imply a higher duty cycle of nuclear "delta" -bursts. It 
thus appears that, at least to some degree, the mechanisms that regulate the (past) star formation 
efficiency in the vicinity of the galaxy nucleus - as traced by the central surface brightness of its 
stellar disk - also play a significant role for the amount of gas that is funneled into the very nucleus 
and converted into stars. Further support for this notion comes from the fact that the molecular 
gas content within the central kpc of late-type spirals also appears to be connected to the disk 
surface brightness in the sense that low surface brightness disks in general have little or no central 
molecular gas (Fig. 7 of B6ker, Lisenfeld & Schinnerer 2003b). 

4.3. Hubble-type 

Figure 11a compares nuclear cluster absolute magnitudes against the Hubble type of their host 
galaxies. The very weak correlation apparent here, with a Spearman coefficient s = —0.26, is only 
marginally significant (at the 2-a level, and less if the cluster in NGC 7418, with Mp = —16.33, 
is excluded). It seems in any case to be the result of the strong link between NC and total galaxy 
luminosities, together with a trend towards slightly fainter galaxy magnitudes, on average, for the 
later Hubble types. Figure lib therefore plots the ratio of cluster-to-disk luminosity vs. Hubble 
type. The correlation coefficient in this case is only s = 0.13, which is not significant at even the 
l-cr level. 

Thus, nuclear cluster luminosity is essentially independent of Hubble type. This should not 
come as a surprise given that the Hubble type is not a strong discriminant for the latest-type 
spirals, as discussed in more detail by Boker, Stanek, & van der Marel (2003). The main reason 
is that the NC is quite easily mistaken for a compact stellar bulge in seeing-limited images, and 
hence a ground-based morphological classification based at least in part on the bulge-to-disk ratio 
is inaccurate at best. 

4.4. Nuclear Clusters and Stellar Bars 

Using the morphological classifications of our galaxy sample, both from the RC3 (listed in 
Table 1 of Paper I) and the LEDA database, we have investigated for our complete Paper I sample 
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whether nuclear star clusters are more likely to be found in (or whether they preferentially avoid) 
galaxies of a certain bar class. Table 3 summarizes the statistics for both subsamples with (59 
galaxies) and without (18 galaxies) a NC. The fraction of galaxies with a NC is 76%, 88%, and 67% 
for the RC3 classifications A, AB, and B, respectively. The LEDA database only lists the types 
"unbarred" and "barred"; for these, the corresponding fractions are 88% and 73%. The apparently 
smaller frequency of NCs in barred galaxies is probably not significant, given the small number 
statistics in some categories, and the rather high number of uncertain classifications (17 galaxies in 
our sample are classified either as uncertain or peculiar in the RC3). More reliable classifications of 
the morphology of late-type spirals will require a more accurate isophotal analysis or even kinematic 
information on their central disks. 

This rough analysis suggests that the galaxy bar class is not a significant factor for the likeli- 
hood of hosting a NC. In particular, the late-type galaxies in our sample seem to have no especial 
difficulty in supporting bars and nuclear clusters simultaneously. This might have seemed some- 
what paradoxical in the light of some early numerical simulations of barred galaxies, which found 
that fairly modest central mass concentrations could drive dynamical instabilities that destroyed 
the bar on rather short timescales (e.g., Friedli 1994; Norman, Sellwood, & Hasan 1996). However, 
more recent work (Shen & Sellwood 2003) suggests that even an NC as massive as a few percent 
of the total disk mass of a galaxy (i.e., with C/D > 0.02, already much higher than the typical 
values in Fig. 10b) may not be able to destroy a bar completely within a Hubble time. It is then 
somewhat unclear what role bars and NCs together might play in the formation of the exponential 
"pseudo" -bulges found in many late-type spirals (cf. Carollo et al. 2001). 



5. Summary 

We have presented a structural analysis of 39 nuclear star clusters in late-type spiral galaxies. 
After correction for the instrumental point spread function, we measure a median effective radius 
for the sample of = 3.5 pc, with 50% of the sample falling between 2.4 pc < r e < 5.0 pc. This 
narrow size distribution is statistically indistinguishable from that of Galactic globular clusters, 
even though the nuclear clusters are on average 4 magnitudes brighter than the old globulars. The 
compactness of NCs clearly separates them from larger ellipsoidal stellar systems such as spiral 
bulges or elliptical galaxies. 

In addition, we have compared the luminosity of NCs with various properties of their host 
galaxies. We find that NC luminosity is correlated with galaxy luminosity which extends a pre- 
vious result obtained for earlier-type spirals. Because there is a strong correlation between total 
luminosity and central disk surface brightness for the galaxies in our sample, NC luminosity also 
correlates with local disk brightness. It is not clear which of these trends is physically more fun- 
damental. Either way, other apparent correlations such as between cluster luminosity and galaxy 
size can ultimately be explained by the underlying link to galaxy luminosity and/or central disk 
surface brightness. 
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We find no systematic trends with the Hubble type of the galaxy, confirming that the latter 
is not a strong discriminant for the latest-type spirals. We also find no evidence for a connection 
between the presence/absence of NCs and the presence/absence of large-scale stellar bars in the 
host galaxy. 

Without additional spectral information, it is not possible to put strong constraints on the 
formation history of NCs. We suggest that the most likely scenario for NC formation consists of 
multiple nuclear starbursts, possibly caused by gas infall onto a "seed" cluster. However, a more 
detailed discussion of this hypothesis has to await the analysis of our spectroscopic data which we 
will present in a future paper. 

We are grateful to S. Larsen for help with the implementation and use of his ISHAPE software 
package. Support for this work was provided by the National Aeronautics and Space Administration 
(NASA) through grants for project no. GO-8599 awarded by the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract no. NAS5-26555. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with NASA. It has also benefited 
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Table 1. Summary of Cluster Properties 
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NGC4701 


11.0 


15.91 


16.80 


-13.46 


KinglOO 


0.046+0.001 


2.45 


5.44 


NGC4775 


22.4 


16.82 


18.04 


-13.78 


KinglOO 


0.039+0.004 


4.24 


5.09 


NGC5068 


8.7 


18.21 


17.54 


-12.34 


King30 


0.177+0.064 


7.47 


4.03 


NGC5584 


24.2 


17.68 


21.91 


-10.08 


ngf 1 








NGC5585 


8.2 


17.92 


18.24 


-11.35 


KinglOO 


0.082+0.006 


3.26 


4.35 


NGC5668 


23.8 


17.80 


18.85 


-13.10 


KinglOO 


0.038+0.004 


4.38 


4.79 


NGC5669 


21.2 


18.45 


21.65 


-10.03 


ngf 1 








NGC5774 


23.5 


18.78 


21.96 


-9.98 


ngf 1 








NGC5964 


22.2 


18.43 


19.21 


-12.63 


KinglOO 


0.046+0.003 


4.95 


4.50 


NGC6509 


27.5 


17.39 


19.48 


-13.09 


KinglOO 


0.021+0.005 


2.80 


5.18 
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Table 1 — Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Galaxy 


Distance 


Disk fj° 




Mf 


Best Model 


r e 


r e 


log(/e) 




[Mpc] 


[mag/arcsec 2 ] 


[mag] 


[mag] 




[arcsec] 


[PC] 


[L /pc 2 ] 


NGC7418 


18.4 


14.64 


15.02 


-16.34 


ngf 2 








NGC7424 


10.9 


17.23 


18.79 


-11.42 


ngf 2 








NGC7689 


24.9 


16.49 


18.22 


-13.78 


ngf 2 








NGC7793 


3.3 


17.59 


13.99 


-13.64 


KinglOO 


0.484±0.006 


7.74 


4.51 


UGC3574 


23.4 


18.35 


19.97 


-11.98 


KinglOO 


0.068±0.008 


7.71 


3.85 


UGC3826 


27.8 


19.11 


21.59 


-10.76 


ngf 1 








UGC4499 


12.5 


19.83 


21.91 


-8.65 


ngf 1 








UGC4988 


24.2 


19.23 


20.76 


-11.20 


KinglOO 


0.030±0.006 


3.52 


4.22 


UGC5015 


25.7 


19.14 


20.71 


-11.37 


King30 


0.060±0.017 


7.48 


3.64 


UGC6931 


20.7 


20.15 


21.87 


-9.75 


ngf 1 








UGC8516 


16.5 


18.15 


20.16 


-10.99 


ngf 2 








UGC12732 


12.4 


20.14 


19.35 


-11.29 


KinglOO 


0.065±0.003 


3.91 


4.17 



Note. — (2): galaxy distance from Paper I. (3): observed (i.e. not inclination-corrected) peak surface brightness 
of the galaxy disk (underlying the NC). (4): apparent /-band magnitude of NC. (5): absolute /-band magnitude 
of NC, corrected for foreground reddening. (6): name of best-fitting ISHAPE model, "ngf" means that no good 
fit was possible, either because of low signal-to-noise data (denoted with superscript 1), or a complex morphology 
(superscript 2). (7) angular effective radius and uncertainty as discussed in § 2. (8) effective radius in parscc. (9) 
logarithm of effective I-band intensity of NC, defined as log(/ c ) = 0.4(4.08 — Mf) — log(2-7rr 2 ) for r e in pc. 
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Table 2. Fit coefficients and statistical significance 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


X 


y 


Figure 


a 


b 


SROCC 


Sign. Level 


m% 


m^ 7 


6a 


6.8±1.5 


0.99±0.11 


0.757 


5.762 


Mf 


Mf 


6b 


2.6±2.3 


0.78±0.13 


0.654 


4.978 


M FIR 


Mf 


7 


-2.0±3.7 


0.53±0.20 


0.420 


2.554 




Mf 


9a 


-29.2±1.6 


0.60±0.09 


0.679 


5.173 




Mf 


9b 


-25.2±0.6 


0.75±0.11 


0.627 


4.778 



Note. — Cols. (1-3): quantities to be correlated using the expression y = 
a + b X x, and Figure in which they are plotted against each other. Cols. (4) 
& (5): coefficients a and b of best linear least-squares fit, together with their 
statistical uncertainties (lcr). Cols. (6) & (7): Spearman rank-order correlation 
coefficient (SROCC) and significance (in a) of correlation. 



Table 3. Impact of bar class on nuclear cluster frequency 



(1) 


(2) 


(3) 


(4) 


Database 


Class 


with NC 


without NC 




unbarred (A) 


76% (16) 


24% (5) 


RC3 


mixed (AB) 


88% (23) 


12% (3) 




barred (B) 


67% (20) 


33% (10) 




total 


77% (59) 


23% (18) 




unbarred 


88% (15) 


12% (2) 


LEDA 


barred 


73% (44) 


27% (16) 




total 


77% (59) 


23% (18) 



Note. — Sample fraction and total number of galaxies 
with and without a clearly identified nuclear cluster, sorted 
by bar classification. For comparison, we have listed the 
statistics for both the RC3 and LEDA databases. 
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Fig. 1. — Some typical examples of ISHAPE fits. The three panels of each row show (from left 
to right) the data, the best ISHAPE fit (analytic model convolved with TinyTim PSF), and the 
residuals (data — fit). The field size of each panel corresponds to the aperture radius R u (Paper I) 
over which the NC dominates the emission. From top to bottom, the objects are NGC 2805, 
NGC 3445, NGC 3906, and NGC 4618. The last was excluded from the analysis because of its 
complex morphology which prevented a good fit. 
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Fig. 2. — ISHAPE fit results for NGC 2805. Shown are the x 2 values of each fit as a function of 
effective radius for all analytic models explored. These values have been rescaled such that the 
minimum x 2 matches the number of degrees of freedom (see text). The inlay panel shows that only 
the MOFFAT15 model with r e values of 0.056" ± 0.006" falls below the 3a confidence limit for 5 
free parameters (Ax 2 = 18.2, indicated by the dashed horizontal line). Within the MOFFAT15 
model, the range of acceptable vaues for r e is marked by the solid horizontal line which denotes the 
3cr confidence limit for a single free parameter (A% 2 = 9). 
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Fig. 3. — Comparison of NC effective radius and absolute magnitude. The value of r c is not a 
strong function of cluster luminosity: clusters of similar size span a wide range in magnitude. Two 
clusters with large r e values have been omitted from the plot in order to avoid crowding of the 
data points, they are the NCs in NGC450 (r c = 24.3 pc, Aff = -11.96) and NGC4625 (r c = 28.8 pc, 
Mf= -10.63). 
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Fig. 4. — Left: distribution of total absolute I-band magnitude for all NCs of Paper I (solid line), 
compared to that of the Milky Way GCs (dotted line). Right: distribution of effective radii for 
all NCs with good ISHAPE fits (solid line), compared to the Milky Way GCs (dotted line). All 
histograms have been normalized to a peak value of 1. 
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Fig. 5. — Comparison of NCs with other spheroidal populations in the r e vs. I e plane. The data are 
from the following sources: Galactic GCs: Harris (1996), Ellipticals, spiral bulges, and dE galaxies: 
Burstein et al. (1997), dE nuclei: Geha, Guhathakurta, & van der Marel (2002). All photometry 
has been converted to I-band, assuming the following colors: B— 1=2.0 for E, SO, and dE galaxies 
from the Burstein et al. (1997) database, V— 1=0.8 for the dE nuclei of Geha et al. (2002), and 
V— 1=0.95 for those GCs which have no value for Mj listed in the Harris (1996) catalog (see § 3.1). 
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Fig. 6. — /-band mass-to-light ratio (top) and absolute /-band magnitude for a stellar cluster of 
initial mass 7 x 10 5 M Q , formed in a single burst of star formation with the IMF of Kroupa et 
al. (1993) between lower- and upper-mass cut-offs of 0.1 M and 120 M & , as derived from the 
PEGASE code of Fioc & Rocca-Volmerange (1997). The predictions include the effects of mass 
loss due to stellar winds and supernovae, but not the effects of dunamycal evolution. The four lines 
of each panel correspond to different metallicities as indicated. 
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Fig. 7. — Comparison between NC luminosity and total luminosity of the host galaxy, a) apparent 
I-band magnitude of NC vs. apparent total B-band magnitude of the host galaxy (from the LEDA 
database), b) same as a), but for absolute magnitudes (calculated using the distances listed in 
Paper I). The median error bar is indicated in the upper left corner of each panel, excluding (for 
panel b) any (unknown) distance errors. Both panels show a statistically significant correlation, 
which is indicated by the the straight line in each panel. The lines mark the best least-squares 
linear fit to the data (fit coefficients and significance levels are given in Table 2. 
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Fig. 8. — d) I-band magnitude of NC vs. IRAS FIR magnitude of the host galaxy (from the LEDA 
database). Over the luminosity range of our galaxy sample, only a weak correlation exists. The 
lines mark the best least-squares linear fit to the data. The fit coefficients and significance levels 
are given in Table 2. 
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Fig. 9. — a) I-band magnitude of nuclear star cluster, Mp, and b) total blue galaxy magnitude 
Mg, as a function of central disk surface brightness n®. All three quantities are correlated with 
each other. 
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Fig. 10. — a) I-band magnitude of nuclear star cluster, Mp, and b) cluster-to-disk luminosity ratio 
C/D as a function of physical isophotal diameter of the galaxy disk. Any apparent trend can be 
explained by the strong correlation between cluster luminosity and galaxy luminosity (see text). 
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Fig. 11. — a) I-band magnitude of nuclear star cluster, Mf, and b) cluster-to-disk luminosity ratio 
C/D as a function of galaxy Hubble type. No trends are apparent. 



